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Abstract 

In this paper, we use canonical transformations to cohectively analytically continue 
the singularities of the simultaneous binary collision solutions for the collinear four- 
body problem in both the decoupled case and the coupled case. And more importantly, 
we describe the relationship between the decoupled solutions and the coupled solutions. 

1 Introduction 

The question of the regularization for a simultaneous binary collision (SBC) solution is not 
completely understood although many results about it have been obtained. In [1] , Saari 
showed that a SBC solution can be analytic continued by rescaling the time s = ta and a 
majorant method argument. In [3], Simo and Lacomba gave a different approach and also 
they showed that "simultaneous binary collisions in the classical n-body problem are 
block-regularizable" . In p], Elbialy proved that "collision and ejection orbits can be collec- 
tively analytically continued, i.e. each collision-ejection orbit can be written as a convergent 
power series in ts, with coefficients that depend real analytically on the initial conditions". 
In [1], Martinez and Simo also discussed the block regularization and the result is this "reg- 
ularization is differentiable but the map passing from initial to final conditions is exactly 
(^8/3" _ Punosevac and Wang constructed coordinate transforms that removed the sin- 

gularities of simultaneous binary collisions in a pair of decoupled Kepler problems and in a 
restricted collinear four-body problem. 

In this paper, we use similar canonical transformations and time transform s = s{t), which 
were used by Siegel and Moser([6]) to do the analytic branch regularization at the singular 
point for a binary collision. To the best of our knowledge, this is the first time canonical 
transformations are introduced for SBC. In Section 3, we use the extended analytic im- 
plicit function theorem to show the branch regularization for the decoupled case with special 
masses. And also, we find the special constant C, which is crucial for the decoupled case. 
In Section 4, we use the majorant method argument to prove the branch regularization for 
the coupled case with special masses. In Section 6, the coupled case with general masses are 
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discussed and it follows the same argument as Section 4. In Section 5, we calculate the Taylor 
series solutions for both the decoupled case and the coupled case. By comparing them, lots 
of connections are given. 



Express the equations of motion for the coUinear four-body problem in Hamiltonian form, 
assuming that the center of mass Pq rests at origin. We denote the coordinates of the points 

Pi, P2, P3, P/i by qi, q2, qs, q^, and denote the momenta by pi, p2, p^, respectively, where 
Pk = fTT-kQk {k=l, 2, 3, 4). Without loss of generality, we assume qi < q2 < < on the 
line, pi > 0, p3 > 0, p2 < and p4 < 0. For the coUinear four-body problem, we have the 
Hamiltonian system: 

, E^T-U{k^l,2,?,,A) 

Pk = -J^qu 



rriimj 



l<j<i<4 I I 

Also, assume the center of mass and the total momentum both are o. i.e. 

miqi + m2q2 + ^3^3 + 'm4q4 = 0, 
miQi + 171242 + m^qz + rn^qi = 0, 

or 

Pi +P2+P3+P4 = 0. 

The first canonical transform (Section 4) is defined as foUowings: 

Xi^q2- qi, 
^2 = ?4 - ?3, 

^3 = ^(^2 - qi) + ^(^4 - qs) +q3- q2] 

5^1 = -Pl + ^(Pl+P2), 

^2 =P4+ ^(Pl +P2), 
^3 = -Pl -P2; 

Followed by a second canonical transform: 

^l^-XiYl i2^-X2Yi, 6 = ^3, ^1 = ^, V2^^. 773 = 13 

-'1 -'2 



First, we work on the decoupled case: ^3 = 00, rjs — h — 0. In proving the existence of 
analytic solutions of the new variables for this special case, a constant C = lim^^o ^^l^+llv^ 
appears from our investigation. This constant plays an important role in the existence of 
analytic sohitions and also it gives us a onc-paramcter family of solutions as s close to the 
singular time. For simplicity, assume the singular time to be 0. The solutions of the decou- 
pled case have several nice properties: 

1. The solutions ^i(.s"; C) = {2(5; —C) and 771 (s; C) = 772(5; —C); 

2. rji, r]2 are odd functions and ,^1 ,,^2 are even functions with respect to s; 

3. The solutions for s < correspond to the behavior of the four bodies before the SBC and 
the solutions for s > correspond to the behavior of the four bodies after the SBC; 

4. Although we can choose different values for the constant C for the solution with negative 
s and for the solution with positive s, the solution will be analytic only when we choose 
the same constant for both negative s and positive s. In other word, the analytic solution is 
unique for fixed common constant C on both negative and positive sides of s; 

In the coupled case, we use method of majorants to regularize the SBC. But when we 
come to the solution, several good properties can be observed: 

1. The decoupled case is exactly the special case when ,^3 = 00, 773 = and h — for the 
coupled case by matching the constant in two different cases ; 

2. In the solution of the coupled case, two constants are founded: one is d = —\C, where C 
is the constant in the decoupled case; the other is cu which is given by the initial conditions 
of ^3, V3 s-iid h^and shows the effect of the coupling terms to the solutions. And also the 
coupling term dcu will start appearing from the term in the power series form of Xi and 

X2 ; 

3. The motion of the coupled case is very similar to the decoupled case: the solutions can be 
separated into the decoupled solutions and some power series of cu only up to the power s^, 
but the exact solution is NOT a linear combination of the decoupled case with some other 

movement ; 

4. Up to the power s^, the solution ^1, ^2) Vit V2 ai'c still symmetric with respect to the 
new constant d: ii{s; d) = ^2{s; —d) and 7/1 (s; d) = 7/2(5; —d); 

5. Since uj is fixed for given initial values, but d will make the solution to be a one-parameter 
family which is similar to the decoupled case, the analytic solution can only happen if we 
choose the same common constant d on both negative and positive sides of s. 

2 Preliminaries 

(a) Simplified Hamiltonian form 

We will use the center of mass integrals to eliminate one pair of variables pk, Qk from these 
8 differential equations, and we will achieve this by taking all the Pk, qt into new variables 
^k, Uk via a suitable canonical transformation. We set 

Pk = W,^, Xk = Wy^ (A; = 1,2,3,4), (1) 

where W{q^y) is a generating function whose Jacobian determinant | Wy^q^ \ is not 0. We 
wish to set up the canonical transformation so that xi becomes the distance between Pi and 



P2, X2 becomes the distance between P3 and P4 and Xs becomes the distance between P2 
and P3, while X4 remains to be the coordinate of P4, i.e. 



Xi = q2-qi, X2 = q4-q3, X3 = qs - q2, X4 = q4 (2) 
and set the generating function as 

W{q, y) = ((72 - qi)yi + (g4 - 93)1/2 + (gs - 92)1/3 + qm- 

Then | Wy^^g^ \= 1, and (1) gives a canonical transformation. The first equation of (1) gives 

Pi = -yi, P2 = yi - VS: PS^VS- y2: P4 = y2 + y4- 

Therefore, 

yi = -Pi, y2 = -pi -P2 -ps, y3 = -pi-P2, y4 = pi +P2 +P3 + P4. (3) 

Note ^4 = , y2 = P4 since Pi + P2 + P3 + P4 = 0. 

And it is easy to see that xi, 3:2, X3, 0:4 are nonnegative but yi and y2 are negative. 

In (2), (3) we have the desired transformation, which we see is linear. Since, in addition, 
it does not depend on t, the new Hamiltonian system is 

Xk = Ey^, yk = -E^^ (A; = 1,2,3,4), (4) 

where E — T — U is regarded as a function of the Xk and y^. Then 



Pk _ 1 r ^1 , (^1 ~ ^3)^ , (ys - y2)^ , (y2 + y4) 
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2 mi. 2 nil "7.2 ma m4 ^ 

fe=i 

^ 17111712 ^ mims ^ miin^ ^ 17121713 ^ m2m4 ^ 771^1714 



q2 - qi qs ~qi q^- qi qs -q2 q^- q2 q^ - qs 

mii7i2 17111713 miin^ 17121113 m2m4 11731174 

Xi X'l + X3 Xi+ X2+ X3 X3 X2 + X3 X2 ' 

We know y4 = and 

= migi + m2g2 + JTisga + m4gf4 
= mi(x4 — X3 — X2 — xi) + m2{x4 — X3 — X2) + m3{x4 — X2) + 1714X4, 

so that 

_ XiiTii + X3{mi + 1712) + X2{mi + m2 + ma) 
mi + m2 + ma + m4 
Therefore, we only have to consider the system 

= Ey^, yk = -E^^ (/c = 1, 2, 3), (5) 

with 

E = T-U, 



u = 



y _ hy]_ ^ jyi - y^f ^ (1/3 - ^ yf,! 

2 mi 1712 1^3 ' 

17111712 mims mim^ m2m^ m2m4 msin^ 

Xi Xi + X3 Xi + X2 + X3 X3 X2 + X3 X2 ' 



Xi Xi -t- X3 Xi -t- X2 -t- X3 U 

(b) Binary Collision 

Assume t = ti is the collision time. 
Then, xi — ^ and X2 ^ simultaneously when t ti. 
The following result is from the work of Belbruno([8]). 

Lemma 2.1 

hm — = a, where a = [ ja, (b) 

t^ti X2 1TL3 + m4 

and 

lim(g2 - gi)(g'2 - qif = 2(mi + 7712), (7) 

t — >ti 

lim(g4 - qs){q4 - q^f = 2{m^ + m^). (8) 
t — ^ti 

Lemma 2.2 xi?/^ and X2I/2 ^^^^ ^^^^ finite when t ti. Furthermore, limt^tj xiyf and 
limt_^ti X2y2 exist, and also 

2 r 2 2(mim2)2 

hm xiy, = hm xip, = , 

t-^ti mi + m2 

1- 2 1- 2 2(m3m4)^ 

hm X2?/2 = hm X2J94 = . 

t^h t^ti ms + m4 



Proof First, we will show both xiyf and X2?/| are finite. 

mim^ mim4 m2m^ 7712^4 msm^ 

XiU = 17111712 + Xi h Xi h Xi h Xi h Xi , 

Xi + X3 Xi + X2 + X3 X3 X2 + X3 X2 

As t ^ ti, we have xi ^ 0, X2 ^ 0, 

and Xi + X3, Xi + X2 + X3, X3, X2 + X3, are all positive finite. 

. . hm XiU = hm \m1m2 + Xi = 17111712 + 0:77737714. 

t^ti t^ti X2 

Note that on the phase space of Hamiltonian system (4), T — f/ = h, where h is the Hamil- 
tonian constant. 
Therefore, when t ^ ti, 

XiT = Xi{U + h) 17111712 + a?773?774, 

that is, 



1 r ^ {yi - ys? ^ {ys - y2f ^ {y2 + yif , ^ 

-xi \ \ \ 777 1 7772 + am^m^. (9) 

2 nil 1^2 iTT'S ^4 



In particular, Xiyf and are finite at collision. Then, by (6), xiyf and X2?/| are finite at 
collision. 

Then, we will use the boundedness of xiyf and X2?/2 Lemma [27T] to show the existence 
of the limits of them. 
Note that from (7), (8), (9): 

2(mi + m2) (10) 



lim Xi( 


'PL 




? 




'mi 


m2 




lim X2{ 






? 


t->ti 




1714 






(l/3 


-y2f 


+ 



2(m3 + nii) (11) 

r rZ/i , {yi-y^f , (1/3-2/2)^ , {y2 + yif. ^, . . 

lim Xi\ 1 1 1 J = 2[mim2 + am^imi) (12) 

mi m2 ma m4 

and 

Pi +P2 +P3 +P4 = 0. (13) 

By the boundedness, we can see xiPi, xip\, xipf, xip\ are all finite when t ^ ti. 

Because, pi+p2 = —ys, and by the Hamiltonian system (4), 

, ^ mims mim4 rn2m^ m2m4 

(xi + X3)2 (a;i + a;2 + xa)^ (0:3)2 (x2 + 0:3)2 

For T < t < ti, since X3 is strictly positive, there exists a positive constant B, such that 
X3 > -B > 0. Therefore, integrate the above identity, 

^3(^1) - ysir) < ^(^1 - ^) ■ {fnims + mirrii + 77121713 + m2mi). 

Since the right hand side of the inequality is finite, y^iti) is bounded above. So pi + P2 is 

finite as t approach ti. 

:. XipI, Pi + P2 are 0(1) as t approach ti. 

And then 



Consider (10): 



lim Xi (pi + P2) =0, lim XiPi {pi + P2) = 0. 

t — >t-i t — >ti 



Pi P2 ^2 _ r /Pi , Pi Pi P2 N2 



2(mi + m2) = lim a;i( ) = lim Xi[ h 

t^ti nil ^2 ^1 ^2 ^2 ^2' 

= lima;iPi(— + —Y + —9 lima;i(pi + P2f — (— + — ) lima;ipi(pi + P2) 

t^ti mi m2 m2 m2 mi m2 

T 2 / \ 2 

= hm XipA 1 ) 

t^h nil ^2 

Therefore, limf^/, xi»? exists and the value is ^("^i^z) _ 

Similarly, by considering (11), we can get the existence of Imit^ti X2p\ and also 

2 r 2 ^{Tn^m^f 

iim X2y2 = hm X2Pa = • I 

t^ti t^ti m3 + m4 



3 Main Results 



In this paper, we showed the regularization of both decoupled case and coupled case. Also 
we explain the connection of the solutions between those two cases. 

In the decoupled case, the following forms give all the possible solutions, which is an one 
parameter set of solutions. 

e° = -1 + ds' - -dh^ - -dh^ + -^d^s^ + ... 
5 25 1125 

n 1 d ^ d^ . 58d^ 7 

77? = -s+ -s^ H H + ... 

'^2 5 10 1125 

n 1 d . cP 5Sd^ 7 

77, = -s H + ... 

'2 2 5 10 1125 

where d is an arbitrary constant and 



T C 1 6-6 

a — = — hm 



4 ^s^o^ivf + C2vl 
In the coupled case, all the solutions are: 

^ ^2' ^5 15 12 ^ 



2880 63 1050 25 



^ ^ 2^ ^5 15 12 ^ 

' 2 ^5 12^ 4o 35 60 ^ 
6775 3 19 - 2 139 72 58 -3, 7 ^, 9, 



' 1008000 700 2100 1125 

= - + (-id - —)s^ + (—5^ + — + — 
'2^5 12^ 4o 35 60 ' 

. 6775 3 19 J 2 139 j2 58 J3, 7 ^, 9, 

( duj^ d^uj d^)s^ + Oisn 

^ 1008000 700 2100 1125 ' ^ ' 

The above results tell us that ^ 

1. In each of the decoupled case and the coupled case, it has one parameter rf, which is 

an arbitrary constant. From the comparison, we can see those two constants are the same. 

Recall the meaning of C in section 3, 0?=— j = — | lim^^o J^^^^- There is only one way 



to make the solution analytic, which is to choose the same d on the side of s < and s > 0. 
But if we choose di on the collision side s < and choose a different constant d2 on the 
ejection side s > 0, it is still a solution of the Hamiltonian system and it is not analytic 
anymore. Consider Xi = —^irjj, we can see the constant d will show up from the power or 

4 

n. 

2. The motion of the decoupled case and the coupled case are very similar. Up to the power 
s^, the coupled case can be considered as a decoupled case adding another motion which is 
related to the initial conditions: h, and rj^; on the other hand, because of the mixed term 
duj, the coupled solution can NOT be considered exactly as the sum of a decoupled solution 
and a special solution which has nothing to do with d] 

3. Up to the power s\ the solution ^i, ^2^, 771, and ^2 are still symmetric with respect to 
the new constant d: ^i{d) = ^2{—d) and rii{d) = ri2{—d); 

4. In the solution of the coupled case, basically there are two constants: d = —\C, where 
C is the constant in the decoupled case; another one uj is given by the initial conditions and 
it shows the effect of the coupling terms to the solutions, and also the coupling term du will 
start appearing from the term in the power series form of Xi and X2; 

5. Since u is fixed for given initial values, but d will make the solution to be a one- 
parameter family which is similar to the decoupled case, and the analytic solution can ONLY 
happen if we choose the same common constant d on both negative and positive sides of s. 



4 Decoupled case with all masses equal to 1 

Replace time variable t by the new independent variable 



+ -^^)dt, {T<t< h) 

Siegel and Moser([6]) have shown that j'' Udt is finite, then si = J''{^^ + ^^)dt is also 
finite. 

Denote ^ by x'^ and ^ by y^.Then the system (5) becomes 

— l^m2 I m-inu ^Vk^ Vk — ~ 7nim2 , m^jru ^^k /c = 1, 2, 3 (18) 

Set F = mi ma I m.^m :f {E — h) = mi,n2 1 mam ^ fT — U — h), where E = T — U = h. Then for the 

XI X2 XI X2 

solution of system (5) on the energy surface E = h, we have 

1 „ „ 1 



~Enr-,^ F^l, ™„™ . E^i 



m\m2 I m3m.4 ^fe ' Vk m\m2 _i_ msnn Vk ' 

Xl X2 Xl ' X2 

Consequently, for the solution of system (5) on the energy surface E = h, (18) can be written 

as 

x', = Fy^, y', = -F,., /c=l,2,3 



with F — "1-1112 "j" m-^m^ (T — U — /l). 

If Xk and Uk are solutions of (5) on the energy surface E — h, F is a constant with respect 
to s because 

3 
k=l 

For simphfication, assume = 0, = oo, h = 0, and mi — m2 — — 777.4 = 1- 
we choose 

F^i±4- (19) 
To solve the Hamiltonian system, 

y'k--F., A; = 1,2 (20) 

with F = j^^^ , we are going to construct a canonical transform. 

x\ X2 

(a) Relationship between x^ and 

First, let's write (20) into explicit forms: 

_ 2yi _ 27/1X1X2 
^ + ^ 2:1 + a;2 

_ 1/1 + 1/i 1 _ . . 

(ir + i;)'^? xi (0:1 + 0:2)'^' ^''^ 
/ _ 27/2 _ 27/23:ia;2 , . 

_ yl + yl ^_ X, 

(ir + i)'^i X2(a;i+a;2)'^- ^'^^ 
Think y^ as a function of Xk{k=l,2). 



Lemma 4.1 // {xi, 0:2, 7/1, 7/2} «5 ^/le solution for the above system (21)-(24), there exists a 
constant C such that 

F F 
yl^ h C, and yl^ C, 

Xi X2 

xiyf - X2yl 
or 6 = . 

Xi + X2 

Proof By (21) and (22): 

dyi^Vi^ F_ 

dxi x[ 2yix\ 



Then separate the variables: 

-F 



J '2yidyi = J 

Therefore, 



„2 



dxi 



yl^— + C, or xivl-Cxi^F (25) 

where C is a constant, which depends on the initial conditions. 
By a similar process, we have 

y^^— + or X2yl-CiX2^F (26) 

where Ci is another constant, which depends on the initial conditions, too. 
Add (25) and (26) together, by using (19): 

yl + yl^ — + — + C + Ci 

Xi X2 

= yl + yl + C + Ci 

Therefore, Ci = — C. Then we can rewrite (26) as 

Z/2 = — -C, or X2yl + Cx2^F. (27) 

X2 



(25)-(27): 



xiyi - X2y2 - C{xi + X2) = 0, 
xiyj - X2yl _ 

Xi + X2 

(b) Canonical Transformation 

As we know, at the collision time ti, |/i — > cxo and y2 00. It would be nice if we can 
make a canonical transformation without singularity at t — ti. 

Prom the two-body problem, we have a transformation as 77^ = — . Similarly, we want to use 
this part to generate our canonical transformation. 

1 1 

yi = — , y2 = — ■ 

Let y — {yi,y2)^, b{v) — ^m^m^'^ ~ {xi,X2)'^- Assume the generating function 

V — V{x,r)). The canonical transformation is given by 

y = Va,{x,rj), ^ = Vr,{x,rj). 

Hence, 

y = V^{x,ri) = b{r]) 



Therefore, 

Then 

So 

In particular, let g{'r)) — 0. 
Since 



V{x,rj) =< b{ri),x> +g{ri) 



we can write down the canonical transformation as 

6 = -xiyl, 6 = -X2yl, vi = —, V2^ — ; 

yi y2 

,2 ,2 1 1 

And the new hamiltonian system is going to be 

^k-F,,, Vk^-F^, A; = 1,2 (28) 

with 

CMvf + vi) _ vl + vi ^29) 



+ 6^2 ^ + # 



(c) Meaning of C 



Since we only know the behavior at time t = ti ot s — si.We may think s = si as the 
initial time for the hamiltonian system (28). 

Without loss of generality, let si = 0. We will consider the following two differential equations 

Vi^-F^„ 772 = -^^2 

with initial conditions: 

r/i(0) = 772(0) = 0, ^,{0) = = -F. 

(25) and (27) can be rewritten in terms of and 77^: 

-6 + c^ivl = -6 - c^2vl = F, 

then 

^ - ^ f - ^ 

-I + C772' -1-C77i' 

Differentiate (29), 



2 I 2 ' 

r)t + r)i 



and similarly, 



^6 



Therefore, 



V[ - (-1 + Cvir ■ -A-^ (30) 

'/l "T '/2 



ry^ = (-1 - Cnir . (31) 
Actually, C is an arbitrary constant. 

When C = 0, it is the most special case. The solution for the equations (30) and (31) with 
initial conditions 77i(0) = ?72(0) = is 

s 

111= 112 = -. 



Then 



= -F = = -F 

-l+Cril -'^-Cvi 



Therefore, xi = X2 and yi = y2- 

So when C = 0, the motions of this two decoupled system are exactly the same. 

When C > 0, that is C = ""'f^Hf^ > 0, hence x^yf > X2yl 
Consider the initial conditions. Assume Xi, X2, yi and y2 are the initial values for the de- 
coupled system. If Xi < X2, which means the distance between the two object Pi and P2 is 
less than the distance between the two object P3 and P4, then yi > y2 because xiyj > x'2y|. 
That is, the initial velocity of Pi is also greater than the initial velocity of P2. Note that for 
each collision system, the force between the objects only depends on the relative distance, 
and then by Newton's second law, the acceleration for Pi or P2 is greater than the acceler- 
ation for P3 or P4. So it is impossible for these two collisions to happen at the same time. 
Contradiction! 

therefore, when C > 0, we can get 

xi > X2 and yi > 7/2- 



By a similar argument, we know when C < 0, we have 



xi < X2 and yi < y2- 
(d)Apply C to the decoupled system 

First, we want to show that the constant C can also be derived directly from the following 
initial value problem: 



^2 



^1 = 



V2 



V2 



6(0)=6(0) = -F, r/i(0) = 7^2(0) = 0. 



Let f{s) = because F = 



# + # + 6^? + 6%'' 

Lemma 4.2 // {6, 6^ Vi: V2} is the solution for the above system with the initial conditions, 
then f{s) is a constant with respect to s. 

Proof 

ds- di,'^'^ di2'^'^ d^^'"^'^ 8^2'"^' 
Uv! + ril) 266r/ir/i(6 - 6) ^ -Uvl + ril) -266r/2r/?(6 " 6) 



-26r?i(6-6) e2vM + vl) . -26r?2(6-6) eiVfa^ + r/j) 

= 0. 

Therefore, /(s) is a constant. | 



Lemma 4.3 Let f(s) = C2, then C2 = C, where C is the constant in Lemma\4^ 



Proof Prom the initial value problem, 

dvi v'l (6^? + 6^1)^ ^Iv^ivf + v^) Hvf + vi) 

Note that F = -f^, then rjf + rjj = -F • (| + |), so 

^6 _ 26^71(6 -6) 

1 2^/71(^1-^2) 



:!2. 

i2 ' il 



6 «i ■ 



Separate the variable, 
integrate both sides. 



where C3 is a constant. 
By the initial condition. 



Therefore, 



that is 



-f^ + 42% -f" 



^"2 ^ , 

3^ = -p • 277ic?r7i 



6 = 7^"' + 



F ^ 



1 _ C2 2 _ 1 



But from the definition of transformation, we have 
then 

C = C2. I 

Therefore, we can use the constant C in the new hamiltonian system. And note that C 
acts as a first integral in the new hamiltonian system, which does not depend on the constant 
F. Rewrite the system: 

2r?ir/i(^ - ^) _ ^1771 
?i — 2 2 — — zor* — ^, 



V2 



By observation, we have 



Vi6 + ?726 = -i^, (32) 

ryie; + r?2e2 = 0, (33) 
r/l^i' + r/^e2' = F', (34) 
r/i^ie?-r/^r;2'e2 = 0, (35) 

-2CF. (36) 



By (32) and (33), 

iViCi + r/26)' = -i" 
Since when s = 0, ?7i = 7^2 = 0, therefore by integrating both sides, 

Vi^i + V2C2 = -Fs 
By the identities -^1 + C^irjf = - C^r^l = F, 



Vi , ']2 



2 + 7^7^ = ^' (37) 



and also the differential equations of rji and 772 will be 



?7^ + ?7i 



Because for C < 0, the equations are 



^2\2 ^2 



= (-1- 1 c hf)- 



Note that the solutions {771,772} of the above two equations are the same as the solutions 
{^2,^1} of (30) and (31) with positive C. 

Without loss of generality, we can assume that C > 0. 



Lemma 4.4 Let {771,772} is the solution for 



r][ = (-1 + Crjir 



2 I 2 ' 



= (-1 - Cr^lf 



r/2 + r^2- 



Define Ni{s) = <^577i(^) an(i N2{s) = C5772(^). T/ien 

tanh~^(iVi) + tan-^(iV2) = s, 

1 - ATf + 1 + ATI ~ 
Proof Consider the fraction between the above two equations: 

77^_ (-l + C77fr Vl 

Separating the variables and integrating both sides: 



1 771 Itanh-^(C^i) _ 1 772 , l tan-HC^2) , ^ 



2C-1 + C77f 2 ci 2CI + C77I ' 2 d 

where D is a constant. 

By the initial condition 771(0) = 772(0) = 0, 

Simplify the above identity of t^i and 772: 

-^J^ + tanh-^(C^7^i) = - tan-i(C^r^2) (38) 

Therefore, combine (38) and (37), 

tanh-^(C^i) +tan-^(Ci772) = + = Ch. (39) 

i + 07^2 1 — 077^ 

Then it is easy to express 772 in terms of 771: 

772 = C~5 tan[C5s - tanh"^(C5r7i)]. 

LetiVi(s) = C^r7i(^), iV2(s) = C-^mi^)- 

Then the equations of 7^1 and 7^2 can be changed to equations of A'"! and N2: 



Ar2 



TV^ = (1 + 7V|)2 



Nl + A^l 
Therefore, by (39), 

tanh-^(A^i) + tan-^(A^2) = s, (40) 

and 



1-Nf 1 + AT^^ 
(e)Existence, Uniqueness and Analytic Properties 

Consider the equations 

Ar( = (1 - 7V2)2 



(41) I 



2 ' 



with the initial conditions 



A^i(O) = A^2(0) = 0. 

Theorem 4.5 The above system has analytic solutions {Ni{s) , N2{s)) as s approaches 0. 
To prove the above theorem, we need to introduce several propositions. 

Proposition 4.6 Assume Ni and N2 satisfy 

No 

tanh-^(iVi) + tan-\N2) = + 



1 + ' 

and {Ni,N2) = (0,0) is a point on the curve given by the above equation, then Ni is a real 
analytic function of N2 in a small neighborhood of N2 — 0. 

Proof: Basically, we will apply the implicit function theorem. 
Compare (40) and (41), we get a equation of A^i and N2: 

tanh-^(Ari) + tan-i(A^2) = ^\r2 + 



1-Nl 1 + Ni 
or 

tanh-^iVO - tan-^(iV2) + + = 



LetG'(A^i, A^2) = -tanh~^(A^i) - tan-^(A^2) + jz^^ + of course there exists a small 

neighborhood V of (0,0), such that G{Ni, N2) is analytic in V with respect to (A^i, A^2)- The 
Taylor series of G{Ni, N2) at (0, 0) is 

G(N,, N2) = y ^I?-7V2n+i ^ y 2n(-l)" 2n+i 
^ ' ^ ^2n + l ^ ^ 2n + l ^ 

n=l n=l 

We know A^^i and N2 satisfy G{Ni, N2) = 0. When applying the implicit theorem, we will see 
that ^f-(0, 0) = 0. So we need to modify the theorem a little bit. 



Proposition 4.7 (Extended Implicit Function Theorem) 

Denote by G'^_^, the second partial derivative |^ by G%^, and the third partial derivative 
S^byG%,If 

G{0,0) = 0, 
G'^^(0,iV2) = 0, 
G';^(0,iV2) = 0, 

G';;^(o,7V2)^o, 

then there exist intervals I — {—5i, 5i) and J — {—52, 62) and a unique function g, such that 

g:J^I, N2^N,=giN2). 

Proof: Differentiate the Taylor scries of G{Ni,N2) with respect to Ni three times, we 
can find some good properties about the partial derivatives of G{Ni, N2) : 

G'^^{0,N2)^0 

G%^{0,N2)^0 
G'^0,N2)=Aj^0 

Because G'^_^ (0, A^2) = 4 > and G%^ is continuous, there exist a rectangular area R: 
\ Ni \< 61, \ N2 \< 62, such that the closure R cV and 

m= min G%(Ni,N2)>0. 

Since G^^(0, A^2) = and G'^^{Ni, N2) is continuous and strictly increasing with respect to 

G"^(7Vi,7V2) >0 for < N, < S, 

and 

Cj^Ti (A^i, N2) < for -Si<Ni< 0. 

By the above result, G"^^(iVi, N2) is strictly increasing with respect to Ni when < iVi < 5i 
and G'j^^{Ni, N2) is strictly decreasing with respect to Ni in (— 0). Note that G^^(0, A^2) = 
0, 

G'^^(iVi,iV2) >0 for 0<N,<5i 



G'^^{Ni, N2) > for -5,<N,<0 

that is 

G'^^{m,N2)>0 for -5^<N,<5i, iVi^O. 
Because ^(0, 0) = 0, G^^(A^i, N2) > when A^i ^ 0,then 

G'(-5i,0)<0, G{6u0)>0. 

By the continuity of G{Ni, N2), there exists < ^2 < ^'2^ such that when | A^2 |< ^2, 

G{-6i,N2)<0, G{6uN2)>0. 

Consider the intervals / = {—61,61) and J = {— 62,62). For any point N2 in J, the func- 
tion G{Ni,N2) is strictly increasing in /, by the intermediate value theorem for continuous 
function, there exists only one Ni & I such that G{Ni,N2) = 0. That means, for any given 
G J, according to G{Ni, N2) = 0, we can always find only one Ni E I corresponds to A^'2. 
By the definition of function, there exist a function g such that 

g:J^I, N2^N,=g{N2). 

Next we need to show g is unique. 
If there exist 

g:J^I, N2^N,= g{N2) 
satisfying G{g{N2), N2) = 0, {N2 G J), then g{N2) = g{N2) {N2 G J). 

Hence, so far we have proven the existence and uniqueness of A^^i as a function of N2 which 
satisfy G(A^i, A/'a) = 0. ■ 

Proof of Proposition 14.61 By Prop. 3.3, the existence and uniqueness are guarantee. 
If we can show A^i is an analytic function of A^2, then we are done. 



Consider 



- tanh-i(iVO - tan"i(iV2) + + = 



1 - A^f 1 + 

Since each function on the left hand side of the equality is analytic close to 0, we can find 
their Taylor expansions for (A^i, in a small interval of (0, 0): 



00 



^2n + l ^ ^ 2n + l ^ 

n=l n=l 



that is 



^2n + l ^ ^ 2n+l ^ 

n=l n=l 



For simplicity, let 



n=2 n=2 



1^ 1^ 3 ^2n + l ^ 

n=2 



and 



n=2 



by ratio test, we can see that hi{Ni) and /i2(-^2) both are analytic in a neighborhood of 
and the radius of convergence is 1. 

In calculus, we know that when r 7^ 0, (1 + xY is analytic for x e (—1, 1) and the Taylor 
series at is 

k=0 

Denote | ^fTT^i""' by miN,) and | ^^^^^^"^ by u^iN,), 

is an analytic function of Ui. Because the composition of two analytic functions is still 
analytic, then [|/ii(A'^i)] 3 is analytic for A^i in a small neighborhood of 0. So does [hi{Ni)]3. 

Similarly, [/i2(-^2)]^ is analytic for A'2 in a small neighborhood of 0. 
Because 

Nf ■ h{N,) ^ Nl ■ h2{N2) 

do the cubic roots both sides. 

By the above argument, both sides are analytic. Let 

r(iVi, N2) = N, ■ [/ii(iVi)]3 - N2 ■ MN^)]"^, 

then r(iVi, N2) is analytic with respect to (A^i, A^2) in a small neighbor hood of (0, 0). 
In order to apply the analytic implicit function theorem, we need to check the conditions: 

r(o,o) = o, 
fir 1 

-(0, 0) = [hi{N,)]-s + N, ■ -[/ii(7Vi)]-i • h[{N^) |^,=o 



,2.1 1 ,2,2 
,3)'+0.-.(-)--0 

= (|)Mo. 



by Cauchy's analytic implicit theorem, there exists ro > 0,and a power series 

00 

N,{N2) = J2^iNi 



1=0 



such that Ni{N2) = J2ilo o-i^l is absolutely convergent for | Ai'2 |< ro and r(A'"i(A^2), N2) = 0. 
That is, Ni is an analytic function of A''2 when | A^2 |< tq. ■ 



Proof of the Theorem 14.51 Since (40) and (41) are true if A^i and N2 satisfy the sys- 
tem. By Prop. 2 and Prop. 3, A^^i is an analytic function of N2 for N2 close to 0. 
By the setting, 

Ni{N2) = ao + aiiVi + a2Nl + ... 

we will show that ag = 0, ai = 1. 
ao = since iVi(O) = 0. 
Because 



^2n + l ^ ^ 2n + l ^ 

n=l n=l 

Substitute Ni by aiA^2 compare the coefficient of on both sides: 

2 , 2 
-a? = - 

3 ^ 3 

then 

ai = 1. 

Comments: since ai = 1, and when s ^ 0, N2 0, 

lim — - = lim — - = 1. 

s^O N2 N2^0 N2 

At the end of this part, we will use the analytic property of A^i with respect to N2 to 
show that both A^i and N2 are analytic functions of s. 
Rewrite the differential equation corresponding to N2- 

N' = (l + N^'^ ' 



2n+l 



or 



'1 + N^f(l- 2 



1 

l + (t)^' 



SetA?"! = X]^o ^11^2 ■, by the claim ^ approaches to 1 as s get close to zero, and when s ^ 0, 
N2 also approaches to 0. so 

^ ^ n=2 



thus, 60 = 0, 61 = 1 and f = 1 + E~ 2 KNr^. 
Set 

N °° 

2 n=l 

where 0(A^2) — Y^^=i c^n-^2 analytic function of N2 in a small neighborhood of 0. 

1 1 



1 1 



2 1 + 

~ 2 



2 

n=0 

which is obviously an analytic function of A^2 in a neighborhood of with radius of conver- 
gence 1. 

Therefore, the right hand side of the above differential equation 

1 



is also analytic with respect to N2 in a small neighborhood of 0. 

By Cauchy's theorem, N2 — {1 -\- N2Y ■ j^r^i N2(0) — has a unique analytic solution 

A^2 = ^2("S) in a small neighborhood of 0. 

And because A^i is an analytic function of N2, then A^i = Ni{s) is also analytic when s is 
close to zero. 

Because r)i{s) = C~^Ni{C^s) for i = 1, 2, then 771 and 772 both are analytic in a neighborhood 
of 0. 

Since ^1 = -i+crp' ^ analytic function of 771 in a neighborhood of and 771 is also analytic in 
a neighborhood of 0, then ^1 is analytic in a neighborhood of 0. And by the same argument, 
^2 is also analytic in a neighborhood of 0. | 

In fact, we can write down the first few terms of the power series solutions of Ni and N2 
for the differential system: 

N[ = (i- Nif ■ ;^ 

K = {l + Nl)^ ' 



with the initial conditions 

7Vi(0) = A^2(0) = 0. 

By an easy calculation, we can get 

Ni(s) = -s - —s^ + —s^ ^^s^ + .. 

^ ^ 2 20 160 36000 

N2(s) = -s + —s^ + —s^ H + .. 

^ 2 20 160 36000 



5 the original system with all masses equal to 1 

By the definition in section 3, the Hamiltonian 

= 1 A^yf , (yi - y^Y , iy^ - y^f y^. 

,m\m2 rriims mim4 m2ms 17121714 17131714, 
+ + + + + — — - h) 

Xi Xi +X3 X1+X2 + X3 X3 X2 + X3 X2 

For simphcity, assume mi = 1712 — — 1714 — 1, then 

F = ±^± iyl ^yl + yl- ym - y2y?,) 

1,111 1 1 1 

-( h — + — + + + + h) 



^ + ^ Xi X2 X3 Xi + X3 X2 + X3 X1 + X2+ X3 

yl + yl yi + y2 1 ^1, 1 , 1 , 1 , ^ 21 . 
= 1 , 1 1 , 1 • ys - 1,1 — + — , — + — , — + — , , — + ^ - z/3 - 1 

— + — — + — — + — X3 Xi + X3 X2 + X3 Xi + X2 + X3 

XI X2 XI X2 XI X2 ±10 ^10 ±1^10 

Right now we will do a canonical transformation such that ■ yz can be absorbed 

XI X2 

into 

x\ X2 

Let ^ ^ 

Yi^yi- 2^/3, ^2 = y2 - -ya, ^3 = ya 

and the generating function 

W{x^, X2, xs, Fi, Y2, Ys) = xi(yi + ^Fs) + X2(Y2 + ^Ys) + xsYs, 



satisfying 



then 



dW , dW 

^ = and ^^X, 



dW 
dW 

""'=972='' 
dW 1 1 



Under the above transformation, the new hamiltonian is 



-[- 



ir + i ir + ^ - 1^1 - 1^2 



-^^3 + h^i 



+ 



+ 



\Xi X^ + hXi + 



+ h- -Y^] - 1 



Let 



A = A{X,,Y:,) 
1 1 1 ^ 



Xz — hXi — hX2 X^ + hXi — hX2 X3 + hX2 — hXi X3 + hXi + hX2 



2 2 ^ ' 2 2 ^ ' 2 ^ 2' 

then 

y2 , y2 1 

X\ X2 Xi X2 

Prom the above Hamiltonian, we can get 6 differential equations: 



X[ = Fy, = 

XL — Fy„ — 



X1+X2 

2Y2X1X2 
X1+X2 
13X1X2 
X1 + X2 



Y'--F -_JI±IL 1 I 1 1 , I 1 , 

= {-F - 1) • ■ ^ + ^^^xi 

1 1 1 



= -i^x. = (-F - 1) ■ ■ ^ + 



X1 + X2 

Follow the idea of section 3(a), we have 

flY. Y' ~ ■ ■ + ^i-L^Xi 



2Fi 

Therefore, by separating the variables, and integrating on the solution surface F 

j 2Y,dY, = j (-^^ + Ax,)dXi 

+ J Ax,dX, + Ci = ^^+ jAx,dX, + Ci 



Where C\ is a constant with respect to Xi . 



By a similar process, 

= + j Ax,dX2 + C2 = ^^+ J Ax,dX2 + C2 

where C2 is a constant with respect to X2. 
(a)New transformation 

By the similar canonical transformation we defined in section 3(b): 

C^^-X,Y,', i^^-X^Yl 6 = ^3, Vi-^. V2^^, m-Ys 

J-l J-2 

therefore, 

ii^2Vlvl ri. 1 2 , 1 



2?l'/l - 2?2'/2 ^ ?3 2'Sl'/l ~ 2?2'/2 ?3 
1 



and the differential equations are 



^2 - --f'a - 9 , t „2\2 + ^2 



(eir/? + i2vi? 

Since we know, when s — > 0, ^3, 773, ^1 and ^2 approach to nonzero constants. In the above 
equations, we can see Mi, M2 , A^i and N2 are 0{s); Gi, G2 are O(s^). 



Introduce a new transformation 



^Ui, ^ = l,2 

s s 2 

= fa + Us, V3 = r)3 + V3, 

and 

s = e~'^, ds = —sdr, 

where ^3 and r/3 are the hmits of ^3 and rj^ ai s = 0. 
Then we can get a differential system about Ui and vf 

dui dvi 1 

-^--F,,+m, _ = + + 

dU2 dV2 1 

-j^ = -F^+U2, ^ = ^^^ + ^^ + 2' 

dr dr 



and 

ds 

dr 

(b)Limits of Ui and Vi at s — 0{i — 1, 2) 



-s. 



Lemma 5.1 



hm Ui — hm U2 — hm Vi = hm V2 — 

S-+0 s-»0 s->0 s-^O 



Proof According to the discussion in section 2, we know 

r/l 1/^ xipf 2(mim2)^ (7713 + 7724) 

hm ^ = hm — 5- = hm — ^ = iim ^ = ■ — . 

s^Q Tjf t^ti Y2 t^ti ^2 ^iPi (^1 + ^2) 2q; (77737774)"' 

Since in our case ttii = 777.2 = 7773 = 7774, 

hm4 = l- 

On the other hand, wc know when t is very close to ti, 7/1 and 7/2 are the same signs. Then 
Iim, positive. Therefore, 

lim^ = l. 

s-^o ni 

By L'Hospital rule, 



hm ^ = hm r)[ = lim ffj"^' + H + hm G^ 



If the limit on the right hand side is finite, then the limit on the left hand side also exists 
and equals to the same value. 

According to section 2, we have 



and 
So 



lim?7i = lim?72 = 0, 



limfi = lim^o = — 1- 



And it is obvious that 



Therefore, 



hm Gi = 0. 



lim — = - . 

s^o s 2' 



and 



Here we can say 



y V2 y V2 y Vl 1 

hm — = hm — • hm — — -. 

s-*Q s s-*o rji s-*o s 2 



lim^;! = lim = 0. 

s-*0 s-»0 



To consider the limit of Ui, we need to go back to Xi and Fj. 
Because hms^o-^i^? = 1, and Yf^ = ?7i = 0{s), then Xi = O(s^). 

Since A is analytic and finite at s = 0, Ax^ is also analytic at s = 0. Consider the integral 
on the interval [0,so], where sq is a small positive number. It is obvious that X[ = 0{s) 
which is bounded on the interval [0, Sq]. And also Ax^ is also bounded because it is analytic 
at s = 0. Hence AxiX[ds — lim£_^o ° Axj^X[ds is bounded by some constant Mq. 
On the other hand, we have 

n' = ^ + y" Ax,dX, + Ci 
X^Y^ -^ = XiJ Ax.dXi + C,Xi 

+ J Ax,dX, + CiX, 

S J S 



that is 



s 

Integrate on [0, sq], and we can see that 



-lim<£l±il 



lim— [lim / Ax,X[ds + C2\ 

s^O S s^O 



Here constant C2 depends on the choice of sq. Then hm^-^o Axj^X[ds+C2 — Jq° Axj^X[ds+ 
C2 is bounded. Note that Xi — 0{s^), therefore, 

hm = 0, or Imvui = 0. 

s-»0 s s-»o 



Prom the equations, we can see 



hm ^ = — hm — = —1, 

^[ s^o r]2 



s 



and 

by L'Hospital's rule. 



hm(Ci + 1) = hm(6 + 1) = 0, 

s— >0 s— >0 



hm — = hm || = —1. 



Therefore, 



\1mu2 — hm — = hm ~^ — • hm — = (—1) -0 = 0. I 



(c) Analytic solutions of and at s = 



So far, we've got a system of 6 differential equations with initial conditions Wj(0) = 
Vi{Qi) = 0. 

dui dvi 1 

du2 _ dv2 „ 1 

^ _ p ^ - _F 

Let s = e~'^, this system can be rewritten as an autonomous system with seven variable Ui, 
Vi and s: 

dui ^ dvi ^ 1 

^ = —^F,, + v^ + -, 

dU2 dV2 I 

-^--F,.+U2, — = F,, + V2 + -, 



and 



dr 



(is 
dr 



For simplification, we may use different notations: 

dak 



dr 



^i=ihi(^i + ^k, {k = l,...,7) 



where a = ((7i, u-jY = {ui, U2, vi, V2, its, "^3, s)^- 

The initial value is 0"^ = 0(A;=1,...,7) and 9?^ are power series in ai, ...,(J7 beginning with 

quadratic terms, and the bki are real constants. 
The seven-by-seven matrix (bki) has the structure 






1 














cu 


1 

























2 


-1 

















-1 


2 







































































-1 



where u ^\h-\rjl^ i. 



Theorem 5.2 The system 



da 

-s— ^ Ba + (p, 
as 



has the initial condition a — (cti, (77)^ = and 



B 






1 














UJ 


1 

















UJ 








2 


-1 

















-1 


2 







































































-1 



where uj — \h 



8% + 



And also ipk(k—l,2...,7) are power series in a\, ...,G'j beginning 

with quadratic terms. 

Then this system has analytic solution a for s close to 0. 



It is easy to see that the eigenvalues of B are —1, —1,0, 0, 1, 1, 3 and B can be dialogized 



as 



-1 00000 

-100000 

3 





1 

1 



So by a linear transformation 



-1 -u 10 

1 1 

-1 1 

r= 10 1 

1 

1 

1 

which satisfies T^^BT — R ,a — T{pi, p^)^ = Tp and {xk) 
can be changed to be 



T ^{(fik), the system 



dpk 

-J— = rkkPk + Xk, 
dr 



(fc = l,...,7) 



where Xk are also power series in pi, ... p-j beginning with quadratic terms. 
Next, we will show that the above differential system 



dpk . f . 

= Jk[P) = TkkPk + Xk, 



{k^l,...,7) 



(1) 



has analytic solutions, where p — (pi, .... py). 

To find the solution, we will carry out substitutions of the special form 
Pk ^ Pk - (pk{pi, P2) (/c = l,...,7) 



(2) 



where (pk are formal series in the first two variables pi, p2 only and begin with quadratic 
terms. If one sets 

3k{p) = jkipi, ■■; Pr) = Xfc + rkk(f>k - (pkpji - (Pkp2f2, {k = 1, 7) 

where on the right p can be expressed as a function of p by means of the substitution inverse 
to (2), then (1) becomes 



dpk 
dr 



rkkPk+jkip), (A; = l,...,7) 



(3) 



with the power series jk beginning again with quadratic terms. We will now determine the 
coefficients of the (pk so that none of the series ji, jr contain product of powers of pi, p2 
alone. In other words, the equations 



jfe(pi,P2,0,0,0,0,0) = 0, (A; = l,...,7) 



(4) 



are to hold identically. 

By (2) the pi, p2 are invertible power series in the two indeterminate variables pi, p2 only, 
and moreover, for ps = 0, py = we have 

Pfe = 0fe(pi,P2), (A; = 3,..., 7). (5) 
Consequently, (4) reduces to the requirement that the equations 



Xk + rkk(f>k - (pkpji - (Pkp2f2 = (A; = 1, 7) 

or 

-rkk4>k + <PkpiriiPi + (f)kp2r22P2 = Xk - <PkpiXi - <Pkp2X2, {k = l,...,7) (6) 

be satisfied identically in pi, p2, where ps,..., p^ are defined by (5). Conversely, from (2), 
(5), (6) we again obtain (4). We now undertake comparison of coefficients in (6). If ap^^ P2 
is a term of 0^ with g(i + g(2 = m > 1, the comparison gives 

{-Tkk + g\r\\ + 92r22)0L = 7 

where 7 is got from a polynomial in the coefficients of the terms in 0i, ... ,^7 of degree less 
than m. Since rn — ri2 — and m > 1, 

-Tkk + giTii + g2r22 = -r-kk -m = l- my^Q, k = 1,2 

For k — 3, 7, rfefe > 0, then of course —Vkk — m 7^ 0. So 

-Tkk + gini + g2r22 = -rkk - m 0, (/c = 1, 7). (7) 

Therefore, induction shows that (4) has exactly one solution in power series 0i, ... , ^7. 
Next, we need to show the convergence of 4>k{k=l,...,7). 

(d) Method of majorants 

Convergence is proved by the method of majorants. 

If 

I I 
are two power series, which need not converge, then g is said to be a majorant of /, symbol- 
ically f -<g,ii 

I Cth.-lm l< K-lm 

for all the coefficients. 
Let 

pi + P2 + ... + p7 = T, Xfc -< _ ^^p , (A; = l,...7). 

Since rn — r22 — —1 and (7) is satisfied, it follows that 

gi + g2 < C2 \ -Tkk + girn + g2r22 \ (fc=l,...7) (8) 

Consequently for the uniquely determined solution ipi, ... ,■^7 of 

ipkpiPi + ipkp2P2 = C2(l + ^kpi + 'tpkp2) 2 _ {k^l,...7) (9) 



Pfc = V'fe(pi,P2) (A; = 3, ...7) 

we have the relation 0^ ^ ■0fe- The reason is following: 
By the previous argument, we have 

-rkk4>k + (pkpirnpi + 0fcp2^22P2 = Xk - 0fcpiXi - 0fep2X2, {k = 1, 7) 

If apf^pf is a term of 0^ with gi + g2 — in > 1, the comparison gives 

(-'^fefe + girn + 5^2^22) a = 7 

then 

I « 1= T 



7 



I {-Tkk + girn + ^'2^22) I 
where 7 comes from the right hand side of (6). 

In the equations of ipk, if /SpiVi^ ^ term of 4>k with g\ + g2 = rn > 1, the comparison gives 



(s-i + 5'2)/3 = C271 

where 71 comes from the right hand side of (9). 
-cir 



Since Xk -< i^cv ~ ■■■'^)' ^^^^ ^^^^ 7i -^1 7 I- Then from (8), 



^ ^ C271 ^ C271 



S'l + 5-2 C2 I -Tkk + + 5'2r22 | 
^ 71 

I -Tkk + S'l^ll + 5'2r22 I 

> 1 ^— r =1 « I .■ 

I -Tkk + S-iril + 5'2r22 I 

By (8), however, ipi — ... — ip^ — ip, and if in addition one sets Xi — X2 — x, it is evidently 
enough to prove the convergence for the solution ip{x) of 

xiIJx = (1 + V'^)T~^^~r~VTT- 
1 - C4[X + ip) 

On the other hand, let il!{x)/x = ip{x), then 

CsX{l + ^f 



(iIj + x^jJ = {l + ijj + xij^)- 



1 — C4X(1 + 1p) 

or 

{i^+xi^J[l - C4X{1+^)] = {1 + + X'^JC3X{1 + 



or 

(■0 + Xlp^) = + 04X1/;^ + Csx{l + Ipy + X^1p^{l + ■0)(C4 + C3 + Csip) (10) 

Let'0 = ctn^"', from (10) we can get the recursion formulas for ak{ k >2) 

ak{l + k) ^ C4ak-i + C4 ^ a^^a^j + 3c3afe_i + 3c3 ^ 

mi+m2=fc— 1 mi+m2=fe— 1 



+C3 ^2 amiamaams + (C4 + C3)(A; - l)ajfc_i + (C4 + 2C3) ^ 

mi+m2+m3=k—l mi+m2=k—l 



+C3 mia^^a^j^ms 

mi+m2+m3=fe— 1 



= [/CC4 + (A; + 2)c3]afc_i + (C4 + 3C3) amiOma + (c4 + 2C3) 

mi+m2=k—l mi+m2=k—l 

+C3 ^ 0^10^207713 + C3 ^ "^la^^a^ja^g (11) 



mi+m2+m3=k—l mi+m2+m3=k—l 

Consider the equation 

*[1 - C4a;(l + *)] = C3a;(l + 



3 



* = C4X-^ + C4X*^ + C3X(1 + (12) 

Let ^ — Yl'^=i ^nX^, from (11) we can get some other recursion formulas for hk{ k > 2): 

bmi bm2 

+ 3c36fe_i + 3c3 ^ 



■'m2 

mi+m2=k—l mi+m2=k—l 



mi+m2+m3=fc— 1 

= (C4 + 3C3)6fc-l + (C4 + 3C3) ^ &mi&m2 + C3 ^ &mi&m2&m3 

m\+m2=k—\ m\+m2+m3=k—l 

Then 



6fe(/c+l) = (A; + l)(c4 + 3c3)6fe-i + (/c + l)(c4 + 3c3) J] 6^16^2 

mi+m2=fe— 1 

^mi ^m2 ^m3 

(13) 

m\+m2+m3=k—\ 

It is easy to see Oi = y and 61 = C3, so ai < 61. 
Assume Oj < 6j for i < k, compare the terms on the right hand sides of (11) and (13): 

{k + l)(c4 + 3c3)6jk_i > kc4 + {k + 2)c3afe_i 
(A; + 1)(C4 + 3C3) 

> (C4 + 3C3) 

mi+m2=k—l mi+m2=k—l 

(C4 + 2C3) Y 



mi+m2=k—l 

(A;+1)C3 Y &mi&m2^'m3 > C3 ^ amiam2am3 

mi+m2+m3=k—l mi+m2+m3=k—l 

+C3 rriiam^am^amg 

mi+m2+m3=k—l 



Therefore, 



bk{k + l) > ak{k + l) 



or 

bk > Ofe- 

By induction, < bk is true for all k. So is a majorant of ip. 

From (12), we can see ^ satisfies a cubic equation. Of course, it has convergent solution. 
Therefore, all the (pk are convergent series. ■ 

By (3) and (4), one obtains for the given differential equation the particular solutions 



Since rn — r22 — —i and s — e""^. 



where di is an arbitrary constant. 
Therefore, the solution for pk is 



4e-'='=-, (A; = 1,2) 
0, (A; = 3,..., 7). 




Pk = ujkipi, A*2) = ujkidie = ujk{dis) {k = 1, 7), 

where Uk are convergent power series in the variables pi and p2 without a constant term, 
and di is an arbitrary constant. 

That is, Pk are convergent series of s in a sufficiently small neighborhood of s = 0. 

Then ak are also convergent power series of s in a sufficiently small neighborhood of s = 

for k = 1,...,7. 



6 Properties of the power series solutions 

^1 = -1 + uis = -1 + ds'^ + {^dw - ld^)s'^+ 

15 5 

17 4 
^2 = -1 + U2S = -1 + {-uj - d)s^ + (--c?2 - —duj - —uj'^)s^+ 

5 15 15 

J16^^3 JLrfa;2 + ILrf2^ + 1(^3)^6 + 0(^8) 

^ 100800 1260 525 25 ^ ^ ' 

s s 3d + jo; 1 2 1 7 1 2\ ^ 

771 = - + Wis = - H ^s'^ + (— H du; a;^)s^+ 

'22 15 40 70 840 ' 



1093 o 1 , o 39 ,o 58 



' 336000 2625 3500 1125 

s s -^u;-3d 3 1 13^ 71 2n 5 

. 6233 ^^^^,_1^^,^_ 58^3),.^ 



144000 250 10500 1125 



^ ^ Iq ^ ^ ^3fil / 611 2 1 7 1 72\ 7 

6^1 60^1 144^1 7^1^ 14400 25 25 ^ 



where 



1, 1-2 1 Iv . 

uj — -h rin + — = - iim A, 

4 8 ^3 4 s^o ' 



^3 

and d is an arbitrary constant. 

If we set (i = (i + ^, then the first four power series solutions can be rewritten as: 
6 = -1 + {d- |).^ + {-\d' + ^/u; - ^u')s'+ 



' 2880 63 1050 25 

6 = -1 + ^)^^ + ^ - ^^^)^^+ 



^ 5 15 

(-- _ 1 J^2 ^ ^^2^ + 1^3)^6 ^ 

^ 2880 63 1050 25 ^ ^ ' 

6775 , 19 ~ 139 5 



' 1008000 700 2100 1125 

' 2^ 5 12^ 40 35 60 ^ 

, 6775 3 19 J 2 139 72 58 J3, 7 ^, 9, 

( o;^ duf (Tuj dns' + O(s^) 

^ 1008000 700 2100 1125 ' ^ ' 

It tells us that each pair Ci = C2, = ^2 when (i = 0or(i=— ^. 
In fact, by the symmetry of the equation, if (^1, ^2, V2, Vs) is a solution of the differential 
system, then (^2, Ci? ^3? V2, ?7i, ^3) is also a solution of the same system. 



Basically, in the solutions, the cocfiicicnts of .^3 and 773 has nothing to do with d, up to 
the power s^. And .^1, ^2, Vi^ V2 has no mixed term du for the first two nonzero terms in the 



power series solutions. 

Compare this solution with the solution for the decoupled case, 

^ 4 80 1600 288000 ' 

^ 4 80 1600 288000 ' 

77? = -S S'^ H S*" s' + ... 

2 20 160 36000 
Q 1 C ^ C'^ 29C'^ J 
^ 2' + 20'' + I60'' + 36000' + - 
where r]^ and ^-i'^ odd functions of s, and rfliC) = rf^{—C)] 

and ^2 are even functions of s, and ^?(C') = C2(~^); 
and F approaches 1 as s approaches 0, which means that F = 1 on the phase space of the 
solutions for the decoupled case. 
If we let C = -Ad, then 

1 1 7 _ 

J4 



-dh^ - 




+ 


7 


5 


25 


1125 




— d^s^ 


+ 


7 





25 


1125 



'^2 5 10 1125 



Therefore, 



n 1 d ^ d^ - 58d^ 7 
^^ = r-r +10' -1125^ +- 
m = V2 - ^'^s' + 0(5^) 

The above results tell us that 
1. In each of the decoupled case and the coupled case, it has one parameter d, which is 
an arbitrary constant. From the comparison, we can see those two constants are the same. 
Recall the meaning of C in section 3, 

T c 1 6-6 

d — = — hm 



4 4 .-0 ^^r]t + 6^2 



2 ' 



2. The motion of the decoupled case and the coupled case arc very similar. Up to the 
power s"^, the coupled case can be considered as a decoupled case adding another motion 



which is related to the initial conditions: /i, ^3 and 773; 

3. Because of the mixed term duo^ the coupled solution can NOT be considered exactly 
as the sum of a decoupled solution and a special solution which has nothing to do with d; 

4. Up to the power sj, the solution ^1, ^2, ?7i, and ^2 are still symmetric with respect to 
the new constant d: ^i{d) — ^2{—d) and r]i{d) — r]2{—d)] 

5. In the solution of the coupled case, basically there are two constants: d = —\C, where 
C is the constant in the decoupled case; another one uj is given by the initial conditions and 
it shows the effect of the coupling terms to the solutions, and also the coupling term dcu will 
start appearing from the term in the power series form of Xi and X2; 

6. Since uj is fixed for given initial values, but d will make the solution to be a one- 
parameter family which is similar to the decoupled^case, and the analytic solution can ONLY 
happen if we choose the same common constant d on both negative and positive sides of s. 



7 Change of Variables for the system with general 
masses 

By the definition in section 3, the Hamiltonian 

F = (T-U -h) 

Xl X2 

^™ + ^«4ml 1112 ^3 

|.mim2 miiTL^ rriim^^ ^2^3 1112111/^ n^m.^ 

Xi Xi + X3 Xi+ X2 + Xs Xs X2 + Xs X2 

™^ + ^^^2mi 2m2^^^ ^2m3 2mJ^^ V2 rris''^^ 

1 1 2 ,^11712 , mim^ mirrii 17121713 17121114 171^1114 

+ {- h - — JUs \ \ \ \ \ h\ 

Zm2 Znis xi xi + X3 xi + X2 + X3 X3 X2 + 0:3 X2 



Let 



mi m4 
yi=yi : 1/3, >2 = 1/2 : 1/3, = Vs 

mi + 1712 + 7714 



Let the generating function is 



W{x,, Yi) = + + X2{Y2 + "2 Y^) + X3F3 

mi + m2 m3 + m4 



So 



Xi = Wy, = xi, X2 = Wy^ = X2, X3 = Wys = — ^ — xi + — ^ — X2 + 0:3. 

mi + 7712 + 1TI4 

And the new hamiltonian is 

^ 1 f mi + m2 2 , ms + m^ 2 

mm2 + m^^2mim2 ^ ^Im^m^ ^ 

^ m2 + ms mi ^4 .-^^2 

2m2m3 2m2(mi+m2) 2m3(m3 + m4'' ^ 
|.mim2 mirris mi 



tl//t2 //ti//t3 //fci//fc4 

^ X3 + ^2^X1 ^X~2 ^ X3 + ^5^X1 + ^5^X2 

^ 7711 +m2 mz+rrii ^ mi+m^ m^+m,4, ^ 



m2m3 

^3 — 

mi+m.2 



^5^X2 X3 ^Xi + ^2^X2 X2 ^ ^ 

1 . mi+m2 2 , ^3 + m4 2n 
+ ^ 2mim2 ^ 2m3m4 ^ 

1 rr rra2 + ^3 mi ^4 ^ 

+ m^il 2m2m3 2m2(mi + m2) 2m3(m3 + m4)^^ 

rriim^ mim^ 

~^X3 + ^2^X1 ^X^ ^ X3 + ^5^X1 + ^5^X2" 

m2m3 m2m4 , _/,]._ 1 

^3 — ; v\i — ; y\2 ^3 — ; Wi -p ; Wo 

mi+m2 m-i+mi ^ mi+m2 m3+m4 

Follow the similar canonical transformation in section 3(b): 

ii^-X,Yl i2^-X,Yl 6 = ^3, m^^, m^^, m-Ys 

J-l J-2 

Xi = -iir^l X2 = -iivl, ^3 = ^3, Yi = -, ¥2 = -, Ys = V3 
and the new Hamiltonian is 

t t (rni+IM.^2 i m3+m4 „2\ 
^ _ ^I'S^V 2mim2 '2 2m3m4 '1^ 

■m^m^iiril + "iim2C2?7i 

^ fr m2+m3 ^ ^1 ^ "-^4 12 ^ ^ 

mziTiAiirjl + mim2i2^l 2m2m3 2m2(mi + m2) 2m3(m3 + m4) ^ 
■miuis mim4 m2m^ 

i: _ 7712 e ^2 I 7714 (T ri2~'~t _ 7^2 C „2 _ rra3 ^ ri^ t _L ttii 2 i 7774 ^2 

?3 7771+7772^1''! ^ 7773+7774 ^2'/2 ?3 mi + m2^1'l m3 + m4^2'/2 ?3 ^^+^3 ?1 '/l ^3_,_^4 ?2'/2 

+ 1-1 

t _|_ 7771 2 _ 7773 2 J 



7771+7772 ^ 'J- 7773 +7774 



The equations corresponding to ^i, ^2, ?7i and 7)2 are 



OC C ^ „2/ m.3m4(mi+m2) /• mim2 (m3+m4) \ 



nC C -^,2^ ;' m3m4(mi+m2) mim2(m3+m4) t \ 



/ ^ _^ ^ mim2^2V2l 2mim2 ^2 + 2m3m4 _^ 
(77737774^7/^ + 777,177126^2)^ 

Consider the hmit of 6 and T/j at s = 0: Similar to the previous argument, we can see 



772 2777?7772 7773 + 7774 , 7773 + 7774,1 

lim = ■ — 5 2- ■ ( )3 

s^O Tjl 777i + 7772 277737771 777i + 7772 



27ni7ni 
hm^ = 

s-»0 7771 + '^2 

hm6= '"^^"'^ 



s-»0 7773 + "^4 ' 

hm77i = hm772 = 0, 

s— >0 s— >0 



lim • '"1+'"^ = hme; = 0, hm ^^^±^ = hm^o = 0, 

s^O S s^O s^O s s^O 

7]i , (7771 +77I2) (7773 +7774)3 

lim — = hm77^ — 



s^O s s^O 277717772(77717772(7773 + 7774)3 + 77737774(7771 + 7772) 



-1 ' 
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T m , (7771 + 7772)3(7773 + 7774) 

hm — = hm 7]^ = j 5- 

s->o s s^o 277737774 [77^17772(7773 + 7774)3 + 77737774(7771 + 7772)3] 

Denote hms_>o ^ — vi and hms_^.o ^ — V2. Do the change of variable 

t I 2mfm| 

_ mi+m2 
1*1 — , U2 



i2 + 


2m'^m^ 






S 


Vl _ 






V2, 


S 




V3 - 


V3 



Vi^ ^^1, V2 

s 

U3 = 6 - 6, ^3 

where 6 a^^d 773 are the limits of 6 and 773 at s = 0. 
Then the new equations become: 

su[ = F^^ - ui, su'2 = F^2 ~ '"2, 



sv[ 



Un 



0: 



-F^2 -V2-V2, 

with the initial conditions at s 

Ui{0) ^ Vi{0) ^ (i = l,2,3) 
Let s — e""^, the above equations can be rewritten as an autonomous system: 



dui 
d7 



-F„^ +ui, 



dvi 
d7 



dr 



du2 
dr 

dv2 
~d^ 



-Fr,^ + U2, 



F^2 +V2+V2, 



and 



ds 
dr 



-s. 



For simphfication, we may use different notations: 
dak 



dr 



The initial value is ak — 0{k—l,...,7) and (fk are power series in ai, ...,(77 beginning with 

quadratic terms, and the bki are real constants. 
The seven-by-seven matrix {bki) has the structure 



B 



'11 


bi2 














bi7 


'21 


622 














b27 








bss 


634 



















644 







































































-1 



where 



^ ^ 2171^1714 y^nii + m2 

■mim2\/ms + + msm4^mi + 1712^ 

2mim|(m3 + 1714) 1 

m3m4 (mi 7712)5 77117772 v^ma + 1714 + 77737774 v^mi + 1712 ' 

2m\ml{mi -\- 7772) 1 

77717772(7773-^7774)! 777177721^7773 7774 77737774 i^777i 7772 ' 



J22 



^33 



1 - 



1 + 



2mim2 v^m3 + m4 



17111712 + m4 + m3m4^mi + m2 ' 
2171^1714 \/m-y + m2 



^34 



mxmisjm-i + 777.4 + m-im^s/mx + m2 ' 



&43 = 



mim2(m3 + m4)i rnxm-i^m-i + 7774 + 77737774^7771 + m2 ' 
—2m\m^{mj, + 77^4) 1 



77737774(7771+7772)5 77717772 v^7773 + 7774 + m^im^^/mx + 7772 ' 
644 



1 + 



2777i7772 v^7773 + 7774 



777 i7772V^7773~+ 7774 + m^m^-^mi + 7772 ' 



&17 = [/i 



777 1 + 7772 + 77l3 + "74 
277727773(7771 + 7772) (7773 + 7774) 

2m\rn\{m^ + 7774 



7^3 + {minis + 77727773 + 777i7774 + 77727774) — 

6 



(777.17772 V^'T73 + 7774 + 77737774^^7777+^772)^ ' 



7771 + 7772 + 7773 + 7774 ^2 , / , , , N ^ i 

527 = Ai - 7; 7 ^ 77 ^ r^?3 + ("^1?773 + 77727773 + 777i7774 + 77727774) — 

27712^3(7711 + 7712) (7773 + 7774) ^3 



2777|777|(777i + 7772) 



(777i7772 V^7773 + 7774 + 77737774 ^7771 + 77*2)^ ' 

To find the eigenvalues of B, we only need to find the eigenvalues for the two different 2 
by 2 matrices: 



hi bi2 

b21 ^22 



Bo 



^33 ^34 
643 644 



By carefully calculation, we find out that the eigenvalues for Bi are 1 and —1; the eigenvalues 
for B2 are 1 and 3. Fortunately, they are exactly the same as the case with equal masses. 
And after some calculation, we can see that B is similar to the same diagonal matrix R: 



R = {rki) 



-1 























-1 























3 







































































1 























1 



Therefore, the previous argument works. And we will also have the analytic properties 
of the solutions of Ui and Vi at a neighborhood of s = 0. 
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